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Nonperturbative calculation of QED processes participated by a strong electromagnetic field, 
especially provided by strong laser facilities at present and in the near future, generally resorts to 
the Furry picture with the usage of analytical solutions of the particle dynamical equation, such 
as the Klein-Gordon equation and Dirac equation. However only for limited field conhgurations 
such as a plane-wave field could the equations be solved analytically. Studies have shown signihcant 
interests in QED processes in a strong held composed of two counter-propagating laser waves, but 
the exact solutions in such a held is out of reach. In this paper, inspired by the observation of the 
structure of the solutions in a plane-wave held, we develop a new method and obtain the analytical 
solution for the Klein-Gordon equation and equivalently the action function of the solution for the 
Dirac equation in this held, under a largest dynamical parameter condition that there exists an 
inertial frame in which the particle free momentum is far larger than the other held dynamical 
parameters. The applicable range of the new solution is demonstrated and its validity is proven 
clearly. The result has the advantage of Lorentz covariance, clear structure and close similarity to 
the solution in a plane-wave held, and thus favors convenient application. 


I. INTRODUCTION 


QED is the successful theory to describe the interac¬ 
tion between particles and photons. Conventionally the 
calculation is carried out in a perturbative manner since 
the interaction is characterized by the fine-structure con¬ 
stant a <C 1. However, in an intense laser field with 
being the four potential, the laser-particle interac¬ 
tion is characterized by the classical nonlinearity param¬ 
eter ^ where e is the electron charge, me 

is the electron mass and <> represents time averaging 
[l[. If $ ^ 1 the interaction is in the nonperturbative 
multi-photon regime. For example, the SLAC E-144 ex¬ 
periment with ^ Ri 0.3 found a electron-positron-pair- 
production rate scaling of i? ^ 0 and in the per¬ 

turbative theory this would be interpreted as the typical 
^ 2 No dependence with iVo = 5 photon absorption. How¬ 
ever, the calculation Q reveals that on average more than 
6 photons are absorbed in the process, and thus demon¬ 
strates the onset of nonperturbative effects on the exper¬ 
iment. It is also found by calculation that if the intensity 
of the laser field in the above experiment is enhanced 
that ^ ^ 1, photon orders up to iV Ri 50 give significant 
contributions to the total rate and thus the process en¬ 
ters the fully nonperturbative regime. Therefore, special 
techniques are required to tackle such strong-field non¬ 
perturbative problems in order to take the effects of the 
laser field properly into account. 

The general approach is to employ the Furry picture 
[^, where the laser field is treated as a classical back¬ 
ground field and the particle state is represented by the 
exact solution of the particle dynamical equation in the 
laser field, so called laser-dressed state. For the ideal sit¬ 
uation of the laser field being a plane wave or a constant- 


crossed field, analytical solutions of the Klein-Gordon 
equation for a meson and the Dirac equation for a fermion 
have been obtained exactly By using such laser- 

dressed state the particle-laser interaction has been taken 
into account to all orders, and the remaining interac¬ 
tion between the laser-dressed particle and the QED vac¬ 
uum is weak, thus allowing calculations in a perturbative 
scheme similar to the conventional QED. Various strong- 
field QED processes in a plane-wave field (or an approx¬ 
imated plane-wave field with the spot radius of the laser 
beam being much larger than the laser central wavelength 
0) have been investigated by this method, including 
laser-assisted bremsstrahlung ll|, multi-photon Comp¬ 
ton scattering 0 , 111111 , electron-positron pair produc¬ 
tion 0, Q, 111 and so on. 

For an arbitrary non-plane-wave field, exact analyti¬ 
cal solutions for the above equations are generally out of 
reach. However, several laser fields of particular inter¬ 
est in study are of this kind, e.g., the field composed of 
two counter-propagating laser waves. The colliding wave 
configuration can provide higher field intensity and the 
electron classical trajectory is distinctly different com¬ 
pared to that in the plane-wave field. Extensive studies 
suggest that this kind of laser field is efficient in pro¬ 
ducing electron-positron pairs and supporting avalanche 
pair production, and thus is ideal for vacuum cascade ob¬ 
servation 161420(1 . Without analytical Volkov solution of 
Dirac equation in this field, the pair production rate has 
been calculated under the important quasi-stationarity 
approximation [l^ : the reaction rate for each pair pro¬ 
duction is calculated with the value of the external field 
fixed at the given time, and finally these rates are time- 
averaged over the total interaction time. This approxi¬ 
mation is justified if ^ ^ 1 since the formation length of 
a QED process in a plane-wave field is ^ times smaller 
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than the laser wavelength P, [^. In 221 a complete 
QED calculation based on Volkov solution in a plane- 
wave field is carried out for the electron-positron pair 
production in the impact of a 50MeV electron with two 
colliding lOkeV X-ray laser beams each of the intensity 
10^°W/cm^. This treatment takes the advantage that 
^ <C 1 though the absolute intensity is high. Although 
it is by this measure not a strong-held problem, novel 
features in the pair production process compared to the 
plane-wave case are identihed. Therefore, it is meaningful 
to investigate analytical solutions for the particle dynam¬ 
ical equations in non-plane-wave helds to calculate more 
accurate strong-held QED reaction rates in these helds, 
especially in the regime with intermediate values of 
The recent breakthrough in solving the Dirac equation 
in pro pagating laser waves is worked out by A. Di Piazza 
I [ 2 ^. Applying the WKB method and looking for 


a solution of the form = ex^{iS{x)/fi\(p(x), the ac¬ 
tion function S{x) is derived hrst by solving the classical 
electron dynamical equation in the background electro¬ 
magnetic held and then the bi-spinor (p(x) is constructed 
via the method of characteristics. In this way the elec¬ 
tron wave functions in the presence of a background elec¬ 
tromagnetic held of a general space-time structure are 
constructed in particular inertial frames where the ini¬ 
tial energy of the electron is the largest dynamical en¬ 
ergy scale. In 2^ it is argued that even in view of the 
future strong laser facilities, such as the Extreme Light 
Infrastructure (ELI) and the Exawatt Center for Ex¬ 
treme Light Studies (XCELS) [i^, it is still necessary to 
employ ultrarelativistic electrons in experimental studies 
on strong-held QED problems in the quantum nonlinear¬ 
ity regime where not only ^ ^ 1 but also the quantum 
nonlinearity parameter y = > 1 with ojh being the 

photon energy in the rest frame of the electron. 


In this paper, we develop a new method and an¬ 
alytically solve the Klein-Gordon equation in a back¬ 
ground electromagnetic held composed of two counter- 
propagating laser waves in a Lorentz covariant manner. 
The method is inspired by the observation that if the 
coefficients of the Fourier expansion of the solution in a 
plane-wave held to different photon modes are written in 
the form of a Bessel function, all parameters can be de¬ 
termined by simple rules, as demonstrated in section 2. 
For the solution in the non-plane-wave held under study, 
the coefficients of the Fourier expansion of the solution to 
different modes of the two kinds of photons are written as 
a multiplication of Bessel functions, and the parameters 
are determined by rules in analogy or as a development 
to the simple rules. The applicable range is obtained by 
examining the validity of approximations used in deriva¬ 
tion steps, and in this way a largest dynamical parameter 
condition is imposed, that there exists an inertial frame 
in which the particle free momentum is far larger than the 
other held dynamical parameters. We can have in mind 
the case of an energetic particle obliquely impacting on 


two counter-propagating optical or X-ray intense laser 
waves in the lab frame, while the calculation can be con¬ 
ducted in arbitrary inertial frames. This is the content of 
section 3. By solving the Klein-Gordon equation, the ac¬ 
tion function S{x) of the solution of the Dirac equation 
is obtained, because if the solution ij; = exp[iS''(a;)/Ii] 
satishes the Klein-Gordon equation 

[{ihdfj, - eAfj,f - ml]ip = 0 , ( 1 ) 

then except for a term proportional to h it results in 

{d^S'+ eA^){df^S'+ eAf^)-ml=0, ( 2 ) 

which is just the equation determining the action func¬ 
tion S{x) for Dirac equation [H, [ 2 ^. 

In section 4 the newly obtained solution is justihed 
by showing that solutions in a plane-wave held can be 
recovered naturally from it. In section 5 the solution 
is simplihed and takes a form in close similarity to the 
solution in a plane-wave held. Finally it is validated by 
substituting the action function into the basic equation 
@ and proving the consistency. 

In the following the natural units c = = 1 is used 

unless claimed otherwise. In the context the electron 
charge e should not be confused with the exponential 
constant in exponential expressions. 


II. RECONSTRUCTION OF THE SOLUTION 
FOR KLEIN-GORDON EQUATION IN A 
CIRCULARLY POLARIZED PLANE WAVE 

We consider the solution for the Klein-Gordon equa¬ 
tion © in a circularly polarized plane wave A = ai cos{k- 
x) + 02 sin(fc ■ x), of which the Fourier expansion to dif¬ 
ferent photon modes takes the form 

OC 


2 p2 

where the dressed momentum is q± = p± + k with 
p± being the particle momentum outside the electromag¬ 
netic held, and Jn is the Bessel function. Our question is 
if we only know the form of the solution is given by Eq. 
m, how to determine the parameters c„, a and /3 so that 
'i/' can satisfy Eq. ([T])? The answer seems to be trivial, 
but the approach to actually obtain these parameters can 
provide useful information. 

Substituting Eq. (IHll into Eq. o, we obtain that 

CXD 

c„ J„(a)e*”^[2nA: • p± + ep± ■ (oi - 

n— — oo 

+ ep± ■ (ai + za 2 )e-''=-"]e“'=-" = 0. (4) 
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It results in the sequence of equations related to different 
photon modes 


c„J„(a)2nfc ■p± + c„_i J„_i(a)e • (oi - m 2 ) 

+ c„+i J„+i(Q;)e*^ep± • (oi + m 2 ) = 0. (5) 

Compared to the identity of the Bessel function 

O' 

7^ —l(o^) “ 1 “ ^n+l(oi)) — (d) 

Zn 

it can be found that in Eq. ([S]) if we let 

Cn = 1, (7) 

and 


Im[e • (oi — m 2 )] = 0 , 


the equation becomes 

—Re[e“*^ep± • (oi — m2)] 
2nk■p± 


( 8 ) 


(fTn — l(o;) “t" (/yi+i(Q;)) — 


and thus 


Re[e • (ai — m 2 )] 

k-p± 


( 10 ) 


From Eq. (jSj) and Eq. m a and /3 can be written in 
the familiar form 


cos /3 = 


p± ■ ai 


sm a = — 


P± ■ 02 


a = — 


ViP± ■ «l)^ + (P± ■ «2)^’ 
e\/(p± • ai )2 + (p± • 02)2 


k- 


P± 


( 11 ) 

( 12 ) 

(13) 


Another solution is got by changing the sign of the 
above quantities, corresponding to /3 = /3 + tt and a = 
—a, but this leads to the same wave function © since 
J„(-a) = (-l)2"e“^J„(a). 

In this section the coefficients of the Fourier expan¬ 
sion of the solution for the Klein-Gordon equation are 
acquired in an empirical way, in which the identity (jS]) 
plays a crucial role. In the following it is shown that the 
experiences and observations gained here can be used to 
tackle nontrivial problems. 


III. CONSTRUCTION OF THE SOLUTION FOR 
KLEIN-GORDON EQUATION IN 
COUNTER-PROPAGATING WAVES 


Our aim is to hnd the solution for the Klein-Gordon 
equation and equivalently the action function for the 
Dirac equation in a non-plane-wave field. Gonsider 

[{p — eA — eA'Y — mg]i/' = 0, (14) 


where A = a[ei cos(fc • a;) -I- e 2 sin(A: • a;)] and A' = 
a'[e'i cos(A:' • x) + 63 sin(fc' • a;)] respectively represent a 
circularly polarized plane wave with the non-plane-wave 
condition k ■ k' ^ 0. As a special case, consider the 
two waves counter propagate with each other. With¬ 
out loss of generality, suppose k = {uj,0,0,kz) and 
k' = {uj', 0 , 0 , fc(,) where kz is positive and is negative, 
and let ei = e( = ( 0 , 1 , 0 , 0 ) and e 2 = e '2 = ( 0 , 0 , 1 , 0 ). 

Based on the observation of the structure of the so¬ 
lution ([3|) in a plane-wave field, we assume the Fourier 
expansion of the solution here takes the form 




E 


E 


C ink-x imk'-x 
nm^ ^ 1 


n— — oo m— — oo 


where 


n — „ J (ly \pinl 3 J I I \ imp' 


(15) 


(16) 


and the formal dressed momentum is written as 




e 2 a '2 


P± =P± 


2k ■ p± 2k' ■ p± 


k'. 


(17) 


The arguments of the Bessel functions are assumed to be 
index dependent, the reason of which will be shown later. 
The quantity p(j_ is not the physical dressed momentum, 
which can only be determined after the acquisition of 
Cnm- The solution m would be obtained by determin¬ 
ing Cnm, am, a'n, /3 and /?'. 

Substituting Eq. (na into Eq. da, a sequence of 
equations related to different photon modes can be de¬ 
rived. Gonsidering the possibility of Cnm to be a function 
of the time-space coordinates, we get 


CnmJniam)Jmian)irQ + TiTi + r2m + rsum) 

+ )Jm{a'n-i)e 

T On+l^mJn+l{,am')Jm{,an-\-\)^ ^b 
T On^m — lJn{,am — l)Jm — l{,an^C ^ d 
T Cn^m+lJniam+l) Jm+l{an)C ^ d =0. (78) 


with 


ro = 


e^a^a'^k ■ k' 


2{k ■ p±){k' ■ p±) 

{p^ — 2eA ■ p — 2eA' ■ p)c 

C'nm 


ri = ±2fc • p± ± 


-t 2e^A • A' 


Cnm 


(p^Cnm)(±p'±), 


ra 


= 2k ■ k', 



(19) 

e2a'2 ^ 2 

: ,, k-k (p^c„m)A:^, 

^ ' P± Cnm 

(20) 

€^CL^ 2 

b - - k ■ k' - iPfj.Cnm)k'''^, 

k * Cnm 

(21) 


(22) 

- *£2), 

(23) 

L - *4)- 

(24) 
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It is found that, the empirical method dzHini) based on 
the identity dH) can be applied. Let 

ro = 0, (25) 

Im(e“*^ 6 ) = Im(e"*'^'d) = 0 . (26) 

Moreover, assuming Cnm is index independent and using 
the relation 

Q! 

(q^) 'An ^TlTTl —— (Oi) -\- — ')tTI 

(27) 


In the derivation for solving Eq. (ESI), we keep h explic¬ 
itly for the benefit of indicating the perturbation orders. 
Presume the solution takes the form 

c„™ = (34) 

where = k — k' and 77 is a constant; /, g and q are 
defined as real functions with the argument k-x, k' -x and 
kd ■ X, respectively. Thus, there is PfiCnm = ifidfj,Cnm = 
Cnm[kfj.f' -f k'^g' (fc^ - k'^)q'] with /', k' and q' being 
the total differential of the corresponding functions with 
respect to the variables fc • a;, k' ■ x and kd ■ x, respec¬ 
tively. Substituting this expression into Eq. EH), it can 
be derived that 


Eq. (IT^ can be written as 

0 — Jni^^m) Jmi^^ri)X\Tl 

H-- [Jn-l(oLm) + d„+i(am)]<4ra((a(j) 

-)- Re(e ^ h^\^Jn — \{^CXrn) Jmis^ri—l) T Jn+l{.^7ri)Jm{.^n+iy\ 
-\- Jnip^m )Jm{a'n)'r2m 

H-['An-l(Q^n) + Jm+l{ot^)\Jn{oLm) 

+ Re(e“*'^'d)[Jm_i(a(,) J„( 

Q^m — 1 ) “1“ {a'n)Jn (o^m+l)] • 

(28) 

If we take the assumption that 

C^mibl ~ and ~ 

Eq. (1281) can be considerably simplified and reduced into 
two equations 

0 = Jn{am)rin 

+ (--h Re(e ^^b))[Jn-l{am) + Jn+l{ctm)], (30) 

0 = Jm{an)r2m 

+ (i:^ + Re(e-^'d))[J™_i(a;) + Jm+M], (31) 


k ■ k'[2f'g' + 2{g' - f')q - 2q‘^ - 2ihq" + QiT 2 g 2 g' 

T 2 ^ 3 /' T 2g'(p3 - Q 2 )] - 2e^aa' cos{kd ■ x) 
T2(f'k + g'k'+ q'kd) ■ p± = 0, (35) 

with gi = 2 g 2 g 3 , g 2 = e^a?/{2k-p±) and ps = e^a'"^/{2k'- 
p±). Eq. (1551) can be simplified by taking f = ±g 2 and 
g' = ±£>3 for corresponding p±, that 

k ■ k'{—2q'^ — 2ihq") — 2e^aa cos{kd ■ x) 

-2{g2k + g^k'±q'kd) ■ p± = (36) 

Therefore, 

± q'kd ■ p± + e^aa' cos{kd ■ x) -I- ^e^(a^ -I- a'^) 

=—{q''^ + ihq")k ■ k'. (37) 

According to Eq. (l37l) ii kd ■ p± = 0, for example in 
a standing wave case = —k^ and a = a' with the 
particle beam shooting perpendicular to the z direction, 
the function 

y{4,) = (38) 

with (j) = kd ■ x/2 satisfies the Mathieu differential equa¬ 
tion 


and thus it is straightforward to use the identity (E]) to 
obtain 


4Re(e-*'^5) 
2ri -I- rsm ’ 
4Re(e-*^d) 
2r2 -I- r^n 


(32) 

(33) 


This result shows explicitly how the arguments of the 
Bessel functions depend on the indices and automatically 
explains the assumption in Eq. EE). 

As a brief summary, the coefficients Cnm, l3, P', am and 
q;(j needed to determine the solution (ITSl) can be obtained 
respectively from Eqs. ([25l [26l [32l and 1331). The remain¬ 
ing task is to solve the Eq. (1251) for the explicit form of 
Cnm and finally check the validation of the assumption 

dMD. 


d(j)' 


2 -I- [ci - 2 c2COs2(/)]?/ = 0, 


with 


2„2 


Cl = - 


2e"^a 


C2 = 


eV 


(39) 


(40) 


A Mathieu equation is also found in solving the Klein- 
Gordon equation in a rotating electric field 27|. The 
formal resemblance is reasonable since a particle locates 
in the vicinity of an antinode of the standing wave expe¬ 
riences a rotating electric field. 

In the following we consider an oblique incidence of the 
particle into the laser field and suppose the condition 


|A| = 


k-k' 
±kd ■ p± 


< 1 . 


(41) 
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Then Eq. (1371) can be solved by the perturbative method 
and the solution is obtained as a sum of terms propor¬ 
tional to different orders of A, that 


q = qo + '^y 


q-n, 


(42) 


n=l 


where the zeroth-order solution takes a simple form 
aa's\n{kd ■ x) [Q 2 k + Qzk') ■ p± 


9o — T" , 

kd ■P± 

Therefore, the solution c. 


"F i. 

kd-p± 

nm reads 


kd ■ X. (43) 


^nm — 

Note that 

g±zp(_''>' 


±i\-g:ik-x-g3k'■x+ ‘ -I- ° ^d'X+O^X)] 


(44) 


^nm — '1^ C , 

(45) 

where the natural units c = h = 1 are used again. 

The condition (I4ip is satisfied in many scenarios. It 
means that there exists at least one inertial frame, here 
for the sake of simple illustration we assume it to be 
the lab frame, in which the particle asymptotic en¬ 
ergy satisfies e S> uj(uj') and we further require that 
e » me^irrieC) with ^ ^ and ^ This is re¬ 

ferred to as the largest dynamical parameter condition in 
the paper similar to that in Ref. 23| . Besides this energy 
requirement, the geometry of the relative movements of 
the particle and the laser waves can become significant 
in special cases, which will be discussed at the end of this 
section. Then it can be calculated that 

ri = ±2fc • p±(l ± ({- -) K. ±2k ■ p±, (46) 

k-p± 

r2 = ±2k'■ p±{lT Ki ±2k'■ p±, (47) 


k' ■ p±‘ 

since the terms -. 

stated earlier, the approximation presented in Eq. 
needs to be checked. It is validated by 


As 


|^m±l 


Wn±l - Oin 


rrie^oj' sin 6p 

e 

vriei' 


£Uj{\ — COsOpY 
me^'uj sin 9p 
euj'Y + cosOpY 


< 1 , 

< I. 


(48) 

(49) 


Therefore, the solution for the Klein-Gordon equation 
(HI is obtained, which reads 


ip = rje 

CO 

X E 




oo oo 

E 

n— — co m— — co 


(50) 


where 


cos /3 = cos /3' = 


P± ■ ei 


sin /3 = sin /?' = — 


V(P± • e-iY + {P± ■ £2)^’ 
P± ■ (-2 


\/{p± ■ eiY + {P± ■ £ 2 )^ 


— 


eaV{p± • eiY + (P± 

k-p±±^k-k' 

ea'\/{p± ■ ei)^ -k {p± - IYY 
k' ■ p± ± ^k ■ k' 


(51) 

(52) 

(53) 

(54) 


and the constant p is acquired in the normalization of 
the wavefunction, that |'0p = 1. 

Finally the applicable range of the solution (IHH) is ad¬ 
dressed. Firstly, let’s scrutinize the steps (HTl) . (HS|) . (ITTl) . 
(H51) and (gSl) which have taken approximations regard¬ 
ing four-vector multiplications. As mentioned previously, 
the validity of these approximations relies not only on the 
largest dynamical parameter condition, but also on the 
geometrical relations among the constituents’ momenta. 
By checking the denominators of the expressions in these 
five equations, that specifically speaking \kd-p±\, \k-p±\, 
and \k' ■p±\, the conclusion is that: to justify these steps 
the azimuth angle 9p of the particle’s free momentum 
direction to the z axis should satisfy 


cos 9p Y 


I / ’ 

OJ + OJ' 


and 


cos0p Y il- 


(55) 


(56) 


Secondly, since the magnitudes of Um and can be 
large since they scale as We^/o; and 1^' respectively, 
the usage of Eqs. (HSl HTl) in Eqs. (15^ 15^ can only 
ensure the neglected part is relatively much smaller than 
the remaining one. To guarantee the absolute magnitude 
of the neglected part to be <C 1, extra constraints on the 
parameters are taken, that 


I ik-k 

\q -|a„ 

k-p± 


^ F’ \q' u ^ l«n < 1- (57) 

k' ■ p± 


These leads to the constraint on the field intensity 




(e^w): 


e'« 


(eV)( 


(58) 


As an example, consider the two waves coming from 
the tunable Ti:Sa lasers with the photon energy oj = 
w' ^ 2eV and the particle being an electron accelerated 
by the laser-plasma wake field to the energy lOGeV [ 2 ^ 
with 9p = 7r/4. Take C ~ 1 which corresponds to 

the laser intensity about 10^®W/cm^. Then it can be 

q'Mi^ - Iq'kT^i 

- oi'nl - 10■^ and W^\a„ 


calculated that |A| - - WM-\ -- 10 




a: 


nil 
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^ 10 and thus it is in the applicable 
range of the solution dSOl). Consider, as another ex¬ 
ample, the same particle condition while the waves are 
IkeV X-ray lasers with ^ ^ ^ 1 corresponding to the 

laser intensity about lO^^W/cm^ which is even higher 
than that available by the present technology, it can 
be found that |A| - 10"'^, W^ W ~ 10"®: 

|am±i - am\ ^ \a'^±i - a'J - lO'^, and \q'^\am ^ 
\^' i^'-p± \ ^n 10“®. Therefore it is also in the applicable 
range of the solution (|5(I)) . 


IV. RECOVERING SOLUTIONS IN A 
PLANE-WAVE FIELD 


By taking A' = 0, it is found that the solution ((5(1)) is 
reduced to the familiar solution of Klein-Gordon equation 
in a plane-wave field as given in Eqs. mB- 

It is more interesting to note that by taking k' in the 
same direction as k and thus denoting k' = rk with r = 
w'/w, the solution dini) becomes 


+ I 1.1 ^ , e^aa'sin(fc,i-x) 




2m/3 imk'-x 


n— — <X) 


m— — oo 


where 


cos /3 = cos j5' = 


P± ■ ei 


sin f3 = sin P' = — 


V(P± • ei)^ + ip± ■ £2)2’ 
P± ■ £2 


a = — 


Vip± ■ + ip± ■ £ 2 )^ 

ea^/{p± -61)^ + {p± ■ 
k-p± 

eaV(P± • £1)^ + iP± ■ £2)^ 


k' ■ p± 


(59) 

(60) 
(61) 

(62) 

(63) 


\ ^ (jj/e in the phase which can be diminished by cal¬ 
culating higher order terms. For the coefficient of each 
photon mode, the approximated arguments lead to an 
error proportional to rrie^fe or rue^'fe. 

Let’s consider the spectral width of the solution. The 
asymptotic formula of Bessel function for n ^ a > 0 
gives 

J„(a)«^(^)", (65) 

V27rn 2n 

Thus the cutoff takes place at n = a and Jn{oi) drops 
sharply as n increases beyond a [ 2 ^. For Jn{a)Jm{o^') 
with 


eav^(p± • ei)2 -P (p± • € 2 )^ 

“ =- Z -’ 

k-p± 

ea' 

^ Uf ' 

k' • p± 


( 66 ) 

(67) 


the cutoff indices are ric = ±a and rric = ±ab These 
turn out to be approximately also the cutoff indices of 
the term Jn{ctm)Jm{pi'n) in tii® solution (l50l) since there 



To measure how well the approximated solution (1501) 
satisfies the Klein-Gordon equation (I14L let’s calculate 
the expectation value of the left-hand side of Fq. (na 
with Dirac bra-ket notation 


=< 'i/'Kp — eA — eA')^ — m\\il} > . (68) 


This is in fact the temporal-spatial integration of the mul¬ 
tiplication of the conjugate of the solution (I5U1) and the 
right-hand side of Eq. (|28l) . It is reasonable to compare 
this residual energy with the characteristic energy of the 
particle, for example in the lab frame or Lorentz in¬ 
variants like k ■ p±. Without going into the details and 
focusing only on the energy scales, we get 




(69) 


Therefore 

with y = k-x and the total field potential At{y) = A{y) + 
A'{ry), and thus the solution of Klein-Gordon equation 
in the two-color plane-wave field is recovered. 

V. VALIDATION AND SIMPLIFICATION OF 
THE SOLUTION 

In the derivation for the solution, we have approxi¬ 
mately solved Eq. (1251) and obtained the arguments of 
the Bessel functions. They are the source of errors to 
the result. In Eq. (1501) . there is an error proportional to 


In the Ti:Sa laser example given at the end of section 2, 
/ 2 11/\2 

~ make the quantity small, one way is to 

reduce the intensity of the laser field to be, e.g., ^ ^ ^ 

O.I corresponding to the intensity about lO^^W/cm^, and 
~ 10“^. In the X-ray example also at the end of 

uju) e ^ 

section 2, there is ~ 10“^. 

It is tempting to write the solution dsni) in a more 
concise form via using the identity 

00 

^ ^ (70) 

— 00 

like what has been done to the solution ((59)) and ((64)) . 
But due to the index-dependent arguments am and a'^ 
of the Bessel functions in the solution (15(11) . not only the 






























7 


direct application of this identity is incorrect, but also 
the sum there can not be decomposed as a product of 
two sums like in Eq. (ESI). This illustrates the complex 
manner of the particle coupling with the non-plane-wave 
field composed of two counter propagating laser waves. 

However, if at the cutoff there is 


- a| - 

1 \ 

^ \aa — ^ 

mUC 

< 1, 

(71) 


e 

uje 


-a'l ^ 

^ \aa — ^ 
e 

^ mUC 

uj'e 

< 1, 

(72) 


which compared to Eq. (ESI) just means the energy scale 
of beje is small, then the following approximation is rea¬ 
sonable 


am ~ a, a„ Ri a . 
Therefore the solution is simplified as 

where the action function reads 


(73) 


(74) 


5 = ± (p± -I- 


-I- e^a'^ 
2 fcd• p± 


kd) ■ x± 


e^aa' sm{kd ■ x) 
kd-P± 


eap±-ei . eap± ■ 62 , 

■ sm(K • X) - - -cos(k • x)\ 


k-p± 


k-p± 


,ea'p±-ei , ea'p±-e 2 , , 

■ sm(A: ■ X) -—- cos(k ■ x)\, 


k' ■ p± 


k' ■ p± 


(75) 


and with <C 1 it is 

{d^S+eA^+eA'^)idf^S+eAf^+eA'f^)-mlil-Oi^O) = «. 

(80) 

Like above, the ratio of the extra term over the charac¬ 
teristic energy of the particle ^ je^ 1. In fact 
there are meaningful strong-held problems with ^ ^ 1 
and 1 which results in <C 1. 

In fact, it is possible to construct other forms of ex¬ 
pressions which can also make the Eq. established 
up to a certain perturbation order, but this paper shows 
how this particular form (1751) of the action function is 
derived step by step and proves its reasonableness. It 
takes the form analogous to that of the plane-wave so¬ 
lution and laser-dressed physical quantities such as the 
dressed momentum and the dressed mass can be directly 
identihed. 

It is worth noticing that neglecting in Eq. (l28l) the 
terms with r^mam and r^na'^, the arguments of the 
Bessel functions would be index independent and the ac¬ 
tion function (ESI) can be obtained. However, if the error 
brought to the Klein-Gordon equation is estimated based 
on Eq. (1^51) . an energy scale of ^ or is found 

which is more severe than found in Eq. (1801) by di¬ 

rect calculation with the explicit expression of the action 
function. 


VI. CONCLUSIONS 


or in an integral form 
S = ±{p±-\ -— kd) ■x± dy 

J —00 


/ 1 2 m/'/' 


2 kd ■ p± 


' —00 

ry 


kd ■p± 




fv I r I 

- / dy'-, -ep± • H(y') - / dy'— - ep± ■ A' 

J -00 k-p± k’-p± 

(76) 

where F{y) = —A ■ A! with y = kd ■ x. In analogy to 
the plane-wave solution, the dressed momentum of the 
particle in this non-plane-wave field can be identified as 


q = p± 


+ e^a'^ 


kd: 


2 kd■p± 

and accordingly the dressed mass is 


(77) 


(78) 


The validity of action ESI can be checked by substi¬ 
tuting it into Eq. ©• Then the left-hand side of the 
equation reads 


{d^S + eA^ + eA'A{d>^S + eA>^ + eH'^) - ' 


-r I 

^3 


{meO 




(79) 


The analytical solution of the Klein-Gordon equation 
and equivalently the action function of a Dirac particle in 
a non-plane-wave electromagnetic field are investigated. 
The method is developed based on the idea that the co¬ 
efficients of the Fourier expansion terms of the solution 
(y')to different photon modes mainly adopt the form of mul¬ 
tiplications of Bessel functions. 

For the field composed of two counter propagating cir¬ 
cularly polarized plane waves, a detailed derivation is il¬ 
lustrated. The coefficients are determined by the rules 
dSSI, (|26l) . (1^ and (|33l as a development or direct anal¬ 
ogy to the rules observed in the plane wave case, see Eqs. 
0 , m and ([TOl . In order to solve Eq. (l25l explicitly 
and justify the assumption (j29p . the largest dynamical 
parameter condition is imposed, that there exists an in¬ 
ertial frame in which the particle free momentum is far 
larger than the other field dynamical parameters. As al¬ 
ready mentioned in section I, this condition is of realistic 
meaning in view of strong-held experimental campaigns 
into quantum nonlinearity regime. The solution for the 
Klein-Gordon equation is obtained analytically, see Eqs. 
(EQHMI. Discussions on its applicable range and exam¬ 
ples can be found at the end of section 3. 

It is clearly shown in section 4 that the solutions for 
the Klein-Gordon equation in the one-color and two-color 
plane-wave helds can be recovered. Considering the cut- 




















off property of the Bessel function, it is found in section 
5 that the solution can be simplified and the action func¬ 
tion takes an integral form (1751) similar to that in the 
plane-wave case. The laser-dressed momentum and mass 
of the particle are identified. The validity of the sim¬ 
plified action is justified by directly calculating the basic 
equation ([2]) that defines the action function. Comparing 
the simplified action d75|l and Eq. (l44|l . it can be found 
that the non-plane-wave feature of the coupling of the 
particle to the two waves is mainly determined by Eq. 

(Eg). 

The solution dsni-o as well as the action Ei is 
Lorentz invariant. It has the advantage of clear struc¬ 
ture and close similarity to the solution in a plane-wave 
case, and thus favors convenient application. With the 
solution of the Klein-Gordon equation, the solution of 
the corresponding Dirac equation in the same electro¬ 
magnetic field can be derived by the method shown in 
paper . It then can provide more exact reaction rates 
of multi-photon pair production process, Compton scat¬ 
tering process and so on in this non-plane-wave field. By 
comparing the results calculated from the plane-wave so¬ 
lution and the non-plane-wave solution, novel features 
particularly related to the non-plane-wave field can be 
identified and used in experimental design and explana¬ 
tion. 

This paper indicates the special convenience of using 
Bessel functions in describing the dressed states of par¬ 
ticles in electromagnetic fields. Extending this method 
to solve problems in other field configurations shall be 
investigated in the future work. 
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